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ABSTRACT 

We report the results of two multi-chord stellar occultations by the dwarf planet (1) 
Ceres that were observed from Brazil on 2010 August 17, and from the USA on 2013 
October 25. Four positive detections were obtained for the 2010 occultation, and nine 
for the 2013 occultation. Elliptical models were adjusted to the observed chords to 
obtain Ceres’ size and shape. Two limb fitting solutions were studied for each event. 
The first one is a nominal solution with an indeterminate polar aspect angle. The 
second one was constrained by the pole coordinates as given by Drummond et al. 
Assuming a Maclaurin spheroid, we determine an equatorial diameter of 972 zb 6 
km and an apparent oblateness of 0.08 zb 0.03 as our best solution. These results are 
compared to all available size and shape determinations for Ceres made so far, and 
shall be confirmed by the NASA’s Dawn space mission. 

Key words: minor planets, asteroids: individual (1, Ceres) - occultations - planets 
and satellites: fundamental parameters 


1 INTRODUCTION 


E-mail: altair08@astro.ufrj.br 


Ceres is the sole example of a dwarf planet in the inner So¬ 
lar System. Far from being mere taxonomic information, this 
suggests the great impact its study can have on the under- 
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standing of planetary formation and evolution of the Solar 
System. Indeed, it was proposed that Ceres’ origin could 
have been as a transneptunian object (McKinnon 20121, 
later scattered to the Main Belt due to the giant planets’ mi¬ 
gration predicted by the ‘Nice Model’ ( Gomes et al.|[2005 |. 
Even if it was formed close to its current location, the dy¬ 
namical history of the Solar System must have left its signa¬ 
tures on Ceres. These could include not only the late heavy 
bombardment features that might exist on its surface, but 
also the makeup of its volatiles, which could have been trans¬ 
ported from the outer region. 

Since the 1970’s it has been speculated that Ceres could 


contain water vapour, which was recently verified ( Kiippers 
et al.|20l4 1. Although the water regime on this object is still 


unknown, some internal structure models suggest the exis¬ 


tence of a water ice - or even a liquid water - layer (Castillo- 
Rogez|20TT (. Yet the very question of whether Ceres under¬ 


went differentiation is open and, on the assumption of an 
affirmative answer, it is natural to ask if it ever had tec¬ 
tonics, what its geological evolution was, and if it is still ac¬ 
tive. Inarguably NASA’s Dawn mission (Russell et al. 20041, 
which is currently orbiting the dwarf planet, will shed light 
on several open issues concerning Ceres. 

Containing approximately one fifth of the whole Main 
Belt’s mass, Ceres is expected to have an equilibrium figure, 
i.e. a Maclaurin or a Jacobi ellipsoid. In fact, direct obser¬ 
vations of Ceres by means of adaptive optics indicate that it 
is an oblate spheroid ( Drummond et ^|2014 l. The precise 
knowledge of its size and shape is of utmost importance, for 
the models of density, internal structure and differentiation. 

The best ground-based technique for determining shape 
and size of a faraway object is the study of its shadow, cast 
by a star during an occultation. Since the 1960’s occulta- 
tions have provided measurements of hundreds of asteroids, 
thanks partially to the fruitful professional-amateur collab¬ 
oration on the field. More recently, this technique has been 
applied to objects of the outer Solar System and has un¬ 
veiled outstanding features of distant bodies, e.g. the ring 
system around the Centaur (10199) Chariklo ( |Braga-Ribas| 
et al.|2014 1. 


The first stellar occultation by Ceres was observed in 
1984 (Millis et al. 19871 and led to the determination of 
its size to the precision of some kilometres, at a time when 
the uncertainties were often ten times larger. The high ap¬ 
parent brightness of Ceres, as compared to most asteroids, 
imposes a somewhat strong constraint on the stars capa¬ 
ble of causing a detectable magnitude drop when occulted. 
For instance, after the 1984 event, to our knowledge, only 
four stellar occultations have been observed ([Dunham et al.| 


2014 1. Two of them had only two chords each, thus not suf¬ 
ficient for providing accurate result^ The two remaining 
events, which occurred on 2010 August 17 and 2013 Octo¬ 
ber 25, are reported in the present work and are the first 
ones that used charged-couple devices (CCD) as recording 
systems. Throughout the paper we shall refer to these events 
as the ‘2010 occultation’ and the ‘2013 occultation’. 

Both events were predicted by Steve Prestoi'0on behalf 


^ These events took place on 1994 August 22 and 2010 October 
30. 

^ Predictions are published at http://asteroidoccultation.com. 



Figure 1. Post-occultation reconstruction of Ceres’ shadow path 
on Earth for the 2010 August 17 event. The big red dot is the geo¬ 
centric closest approach at 22:40:25 UT. The small red ones rep¬ 
resent the centre of the shadow separated by one minute, shadow 
moves from the left to the right. Blue dots are the sites that have 
observed the event. As described in text, UFSC had a negative 
chord. 


of the International Occultation Timing Association, during 
routine prediction of asteroidal occultations of bright stars. 

This work is organized as follows. In Sections 2 and 3 we 
analyse the 2010 and the 2013 occultations, respectively. In 
Section 4 we give the geocentric positions of Ceres derived 
from the occultations. The comparison of our results to those 
in the literature is carried out in Section 5. 


2 THE 2010 OCCULTATION 
2.1 Observations 

As predicted, on 2010 August 17 Ceres occulted the star 
TYC 6833-163-1 (UCAC4 313-111823), which has magni¬ 
tude V — 11.55 and has the ICRS position, based on UCAC4 
catalogue (Zacharias et al.|2013(, to the date of occultation: 


f Q = 17*'18"'29!0085 , . 

\ J = -27°26'38'.'867 ^ ^ 

Observations were carried out in Brazil from five dif¬ 
ferent sites as displayed in Table and Fig. The occul¬ 
tation was detected from four of them. The southern-most 
one (UFSC) had a negative chord. From the positive sites, 
the one named INPE started observing after the event was 
already in progress, due to technical difficulties, thus pro¬ 
viding only the star’s reappearance time; the other three 
recorded the whole phenomenon. 

A remarkable circumstance of this event was the low 
velocity of Ceres: only 3.9 km s”'^ in the plane of the sky. 
Therefore, even exposures of a few seconds would translate 
in a relevant spatial resolution. 

All the observations were made through a sequence of 
images obtained with CCDs. The times of each exposure 
were available on the header of each image with an internal 
accuracy of a few hundredths of a second. CEAMIG had 
only the integer part of the second available, due to the 
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acquisition software used. The fraction of a second could be 
retrieved as described on the Section 3.1 of |Braga- Ribas et| 
(which shall be referred as BR13 hereafter). Cycle 


times (exposure plus read-out) varied from 2 to 52 seconds, 
as can be verified on Table making it a heterogeneous set 
of observations and imposing an error of a few seconds for 
the ingress / egress times of some sites. 


2.2 Light curves 

The flux of the star in the five occultation chords was 
obtained from the FITS images with the platform for 

REDUCTION OF ASTRONOMICAL IMAGES AUTOMATICALLY 

(praia) ( Assafin et al.|20TT |. The light curves were normal¬ 
ized to the flux of the star plus Ceres, as they were merged 
right before and after the occultation. Additionally, they 
were normalized by fitting a polynomial curve (of first or 
second order) outside the flux drop, so that the flux ratio 
was set to 1 outside the occultation. 

The ingress (disappearance) and egress (reappearance) 
instants of the occultation were obtained for each light curve 
by fitting a square-well model convoluted with the Fresnel 
diffraction, the CCD bandwidth, the stellar apparent diam¬ 
eter, and the applied finite exposure time; see|Widemann et| 
[S(] ( [200^ and BR13. The smallest integration time used in 
the positive observations was 1.0 s, which translates to al¬ 
most 3.9 km in the celestial plane. Therefore, the error on the 
time determination of the ingress and egress is largely dom¬ 
inated by the integration times, not by Fresnel diffraction 
or star diameter, which are both on the order of a few hun¬ 
dred meters for this event. The occultation data fit consists 
in minimizing a classical function for each light curve. 


as described in Sicardy et al. (2011) and BR13. The free 
parameter to adjust is the ingress or egress instant, which 


best fittings to the 2010 occultation light curves are shown 
in Fig. and the derived occultation times are listed in 
Table 



Time (seconds after 2010-08-17 00:00:00, UTC) 


Figure 2. The four occultation light curves normalized and ver¬ 
tically shifted by a factor of 0.2 for better visualisation. The solid 
black lines are the best fit of the square-well model to the data. 
Red lines are the square-well model convoluted with the Fresnel 
diffraction, the star diameter, and the applied exposure time. The 
mid-times of the occultations do not coincide due to the propaga¬ 
tion delays of the shadow due to the distinct longitude of the sites. 
Exposures at INPE started after the immersion, as explained in 
the text. 


2.3 Limb fitting methodology 


The methodology used to analyse Ceres’ profile from the 
observations is the same described by [Sicardy et al.| ( |2011| ) 
and BR13. Each combination of site position and recorded 
ingress/egress instant, together with star coordinates and 
Ceres’ ephemeris, corresponds to a point in the plane of the 
sky. The collection of all these points ideally determines the 
apparent limb of Ceres. 

We adopt an elliptic model for the limb profile, resulting 
from the projection of an oblate spheroid onto the sky plane. 
This choice is supported by the work of [Drummond et al.j 
( [2014[ ), by means of direct imaging of Ceres. Hence, we have 
N = 7 chord extremities to adjust the M = 5 parameters 
which define an ellipse: apparent semimajor and semiminor 
axis {a' and b', respectively), position angle P of its semimi¬ 
nor axis and the position (/c, g^) of its centre with respect to 
the occulted star. Of course, the apparent semimajor axis a' 
is equivalent to the equatorial radius Requa of the ellipsoid. 
The coordinates /c and Qc, in kilometres, are calculated us¬ 
ing the JPL#33 Ceres’ ephemeris ( Giorgini et al.|1996 l and 
the occulted star’s position. They are positive toward the lo¬ 
cal celestial east and north, respectively. The position angle 


P is counted positively from the direction of local celestial 
north to celestial east. The apparent oblateness can be de¬ 
fined by e = l — {b'/a'). The best-fitting solution is obtained 
minimizing a reduced Xr function, where we define the num¬ 
ber of degrees of freedom of the problem as Af = N — M. 
All the procedures that allow the determination of the error 
bars of the physical parameters can be found on BR13. 


2.4 Limb fitting solutions 

Two possible solutions were considered for the limb fitting. 
The first, which we call nominal solution, consists of deter¬ 
mining the five parameters that characterize an ellipse from 
the seven observed contacts. As we shortly show, it led to a 
rather large uncertainty on the position angle. Furthermore, 
the nominal solution alone is not capable of returning the 
true oblateness, which can be evaluated through equation 2 
of BR13 provided that the aspect angle is known. 

It is possible to derive the pole position (<(, P) from the 
coordinates of Ceres’ pole {ap,Sp) and its ephemeris (a, (5) 
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Table 1. Circumstances of observation for all observing stations of the 2010 event. 


Site 

Longitude 

Latitude 

Height 

Telescope 

Aperture 

f-ratio 

Exposure 
Cycle time 
Camera 

Result 

Ingress 

Egress 

Observer 

Belo Horizonte 
CEAMIG-REA 

43°59'5T.'1 W 
19°49'4970 S 
825 m 

LX200 

31 cm 
f/10 

5 s 

12 s 

SBIG STIO 

Positive 

22:39:03.9 ± 0.6 s 
22:40:20 ± 5 s 

C. Jacques 

E. Pimentel 

Pico dos Dias 

LNA 

45°34'45'.'l W 
22°32'0377 S 
1864 m 

Zeiss 

60 cm 
f/12.5 

1 s 

2 s 

Andor Ikon 

Positive 

22:37:30.3 ± 0.6 s 
22:41:55.3 ± 0.7 s 

J. I. B. Camargo 
G. B. Rossi 

Sao Jose dos 
Campos 

INPE 

45°51'44'.'0 W 
23°12'3370 S 
620 m 

Cll 

28 cm 
f/6.3 

2 s 

5 s 

SBIG ST7 

Egress only 

Start obs.; 22:39:44 
22:42:03.0 ± 0.2 s 

A. C. Milone 

T. Maldonado 

M. Okada 

Ponta Grossa 
UEPG 

50°05'56'.'0 W 
25‘’05'2272 S 

910 m 

RCX 400 
40 cm 
f/8 

30 s 

52 s 

SBIG STL6E 

Positive 

22:37:17 ± 13 s 
22:39:56 ± 13 s 

M. Emilio 

L. Mehret 

Florianopolis 

UFSC 

48°31'20'.'5 W 
27°36'1273 S 

20 m 

Cll 

28 cm 
f/10 

3 s 

6 s 

SBIG ST7 

No occultation 

Start obs.: 21:49:27 
End obs.: 22:51:21 

W. Schoenell 

A. J. T. Mello 

F. R. Herpich 


cos(C) = — sin((5p) sin(5) — cos{Sp) cos{5) cos(ap — a) 
sin(ap — a) 


tan(P) = 


tan(5p) cos(5) — sin(5) cos(ap — a) 


( 2 ) 

where (" = 0° and C = 90° respectively correspond to a 
pole-on and an equator-on geometry. 

From observations with adaptive optics spanning a 10- 
year period, Drummond et al. (20141 determined the posi¬ 
tion of Ceres’ polar axis within an error of only 3°: 

Qp = (287±3)°, (5p = (-b64±3)°, (3) 


in equatorial J2000 coordinates. Together with Ceres’ 
ephemeris at the moment of the occultation, this corre¬ 
sponds to the polar aspect angle = 86.1°, which is very 
close to an equator-on geometry. Hence, we expect true hg- 
ures to be similar to apparent ones. 

The knowledge of Ceres’ pole not only allows the de¬ 
termination of its polar aspect, it suffices to set its position 
angle. Therefore Eqs. and 1^ may act as a constraint for 
P, and a second solution can be obtained by probing the 
parameter space with the restriction that the position angle 
is confined to the range that follows from Eq. We call this 
the ‘pole-constrained solution’. 


Already mentioned, the parameter with the largest un¬ 
certainty is the position angle: spanning on a 20° interval, 
its determination has a relative precision worse than 10 per 
cent. Clearly, the coordinates of Ceres’ pole (Eq.[^ can im¬ 
pose a strong constraint on the position angle, as the next 
solution shows. 

Finally, the correction to the oblateness due to Ceres’ 
polar aspect angle lies within the 1-sigma error bar and has 
no statistical relevance; hence e = 0.08 ± 0.03. 


2.4-2 Pole-constrained solution 

At the occultation, the coordinates (Eq. of Ceres’ rota¬ 
tional pole correspond to the position angle P = (12 ± 3)°. 
Exploration of the parameter space, restricted to ellipses 
with position angles laying within this range, results in the 
pole-constrained solution. The related physical parameters 
are displayed on the third column of Table in boldface, 
while the best-fitting solution is depicted in Fig[^ 

We notice that the constraint corresponds to the up¬ 
per limit of nominal solution’s 1-sigma error bar for P. On 
the other hand, it selects the smallest values of semimajor 
axis, improving its determination by a factor of about 2. 
Notwithstanding, oblateness’ figures remain the same. 


2 . 4.1 Nominal solution 

With the seven observed contacts it is possible to adjust the 
hve parameters which define an ellipse. For the best-htting 
solution we find Xr,min ~ 0.24, which could be Interpreted as 
a slightly overestimation of the error bars with regard to the 
good quality of the fit. However, inasmuch as the problem 
has only two degrees of freedom, it is far from the statistical 
realm and relatively small Xr,min Eire acceptable. 

The resulting values of equatorial diameter, oblateness, 
position angle and centre coordinates are presented in the 
second column of Table [2] 


3 THE 2013 OCCULTATION 


3.1 Observations 


The event which took place on 2013 October 25 involved 
the star TYC 865-911-1 (UCAC4 496-058191), of magnitude 
V — 10.05. Based on UCAC4 (Zacharias et al. 20131, the 
ICRS position to the date of occultation is: 


f a = ll*'57”52f7641 
\ 5 = -b09°07'49'.'835 
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Table 2. Results of limb fitting to the data of the 2010 and 2013 events. 


Solution 

2010/Nominal 

2010/Pole-constrained 

2013/Nominal 

2013/Pole-constrained 

Equatorial diameter (km) 

982 ± 14 

972 ± 6 

971 ± 7 

971 ± 7 

True oblateness 

0.08 ± 0.03 

0.08 ± 0.03 

0.08 ± 0.04 

0.08 ± 0.04 

Position angle (deg) 

5 ± 10 

12 ± 3 {*) 

22 ± 5 

25 ± 3 (*) 

fc (km) 

97 ± 9 

102 ± 5 

77 ± 6 

78 ± 6 

9 c (km) 

16 ± 15 

21 ± 11 

13 ± 16 

13 ± 16 

■ 

'^r,rmn 

0.24 

0.42 

1.27 

1.27 


Notes. In bold we highlight our best solution. Error bars are at 1-sigma level. The polar diameter (Dp^i) 
can be easily calculated fr om Dppi = Deguaj ^ — e)- (*) Position angle derived from Ceres’ rotational pole 
coordinates determined by [Drummond et al.|(|2014^. 



Figure 3. The best elliptical fit for the occultation chords for 
the event of 2010 using the times from Table ^ and the pole- 
constrained solution. The arrow indicates the direction of mo¬ 
tion, blue lines are the observed chords, the red segments are the 
ingress, egress and mid-occultation error bars at 1-a level. 



Figure 4. Post-occultation reconstruction of Ceres’ shadow path 
on Earth for the 2013 October 25 event at the east coast of USA. 
Upper view of the occultation over the sites that observed the 
event (blue dots). Red points are the centre of the shadow sepa¬ 
rated by 15 seconds. 


The occultation could only be visible in the United States, 
before dawn, as depicted in Fig. 

Nine positive chords were obtained by the variety of 
instruments listed in Table Each station was equipped 
with a video camera with negligible readout times. This was 
of particular importance, since in this event Ceres’ shadow 
speed was 42.6 km s“^, much faster than in the 2010 event. 

Three different timing synchronization procedures were 
adopted among the set of observing stations. At Greenbelt 
and Owings the IPPS signal of a GPS unit was used to 
calibrate time stamps which were inserted at each frame 
of the video. Time extraction is thus straightforward, after 
taking camera delays into account. On the other hand, at 
Brookline the clock would be synchronized by an internet 
server. A lack of connection, however, resulted in spurious 
times. In fact, comparison between the times obtained at this 
station and the others suggests that the former have a delay 
of about 64 s. Therefore, we do not use Brookline’s times 
in the analysis that follows. Finally, at the remaining six 
stations the videos were recorded by camcorders on digital 
tapes. The timing method consisted in the comparison of 
the camcorder internal clock to a IPPS GPS signal, before 
and after the recording of the occultation. Absolute timing 
errors of this procedure are expected to be less than 0.1 s. 

During the event Geres was low in the horizon, with alti¬ 


tudes between 15° (Winchester) and 20° (Hampton). Strong 
scintillation is expected in such a scenario which, combined 
to short integration times and the low magnitude drop of the 
event, resulted in rather noisy light curves and thus larger 
uncertainties in the time of the contacts, as is shown in the 
next section. 


3.2 Light curves 


All videos were converted to FITS images and the photom¬ 


etry of the target was obtained via PRAIA (Assafin et al. 


20111. The light curves were normalized to the flux of a ref¬ 
erence star when available on the field. 

To reduce the noise, the data were binned by groups 
of five images - with the exception of Greenbelt, which was 
averaged in sets of ten. This procedure enlarges the effective 
integration times reported in Table by a factor of 5 (or 
10). As in the 2010 event, an additional normalization by a 
polynomial curve was applied. 

The ingress and egress instants of the occultation were 
obtained by the same procedure explained in Section |2.2| 
Since the typical effective integration time used (0.17 s) 
translates to about 7 km in the celestial plane, and the Fres¬ 
nel scale and star diameter are again of the order of a few 
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Table 3. Circumstances of observation for the observing stations of the 2013 event. 


Site 

Longitude 

Latitude 

Height 

Telescope: 

Aperture 

f-ratio 

Camera 
Cycle time 

Result 

Ingress 

Egress 

Observer 

Hampton 

70°48'59'.'7 W 
42°53'52'.'8 N 

7 m 

12 cm 
f/5 

PC164C-EX2 
0.033 s 

Positive 

09:40:46.9 ± 0.1 s 

09:40:57.26 ± 0.08 s 

T. Blank 

Topsfield 

70°55T6'.'6 W 
42°37'55'.'9 N 
45 m 

12 cm 
f/5 

PC164C-EX2 
0.033 s 

Positive 

09:40:45.4 ± 0.1 s 

09:40:58.0 ± 0.1 s 

T. Blank 

Brookline 

71°08T4'.'5 W 
42°18'27'.'4 N 

109 m 

64 cm 
f/9.6 

Infinity2-1R 
0.015 s 

Positive 

09:41:48.00 ± 0.01 s 
09:42:02.93 ± 0.01 s 

N. Weber 

R. Dantowitz 

Winchester 

78°14'39'.'6 W 
39°16'2T.'5 N 
211 m 

36 cm 
f/5 

PC164C 
0.033 s 

Positive 

09:40:33.26 ± 0.08 s 
09:40:55.86 ± 0.09 s 

J. Brooks 

Greenbelt 

76°52'09'.'4 W 
38°59T2'.T N 
52 m 

12 cm 
f/2.5 

PC164C-EX2 
0.033 s 

Positive 

09:40:33.6 ± 0.1 s 
09:40:56.4 ± 0.1 s 

J. Dunham 

D. Dunham 

Alexandria 

77°02'28'.'3 W 
38°49T9'.'1 N 

8 m 

7 cm 
f/10 

Watecl20N 

0.067 s 

Positive 

09:40:33.3 ± 0.1 s 
09:40:56.1 ± 0.1 s 

P. Maley 

Owings 

76°38'06'.'3 W 
38°41'26'.'5 N 
38 m 

25 cm 
f/3.3 

PC164C 
0.033 s 

Positive 

09:40:34.27 ± 0.05 s 
09:40:56.0 ± 0.2 s 

C. Ellington 

Mechanicsville 

77°23'06'.'7 W 
37°41'26'.'l N 
60 m 

12 cm 
f/2.5 

PC164C-EX2 

0.033 s 

Positive 

09:40:33.0 ± 0.1 s 
09:40:54.8 ± 0.1 s 

D. Dunham 

Varina 

77°19'49'.'3 W 
37°25'58'.'6 N 

19 m 

12 cm 
f/2.5 

PC164C-EX2 

0.033 s 

Positive 

09:40:32.4 ± 0.3 s 
09:40:53.1 ± 0.2 s 

D. Dunham 


hundreds of meters, the theoretical occultation light curves 
are largely dominated by the integration times, as for the 
2010 event. 

The best-fittings to the occultation light curves are 
shown in Fig.[^ and the derived occultation times are listed 
in Table Is] 

3.3 Limb fitting solutions 

Elliptic limb profiles were adjusted to all the availably 
chords by the same procedure described in Section [2.3| This 
yielded Xr,min. ~ 13, suggesting that an elliptic model is not 
satisfactory to the data. Indeed, a quick glance at a plot of 
the observed chords, such as in Fig.[^ shows that the Varina 
chord seems to be somewhat advanced with respect to the 
others. Taking into account that in this station time stamps 
were not inserted on the video frames, it is fairly possible to 
attribute this advance to an eventual problem on the corre¬ 
spondence between camcorder’s and GPS’ times. This could 

® Brookline’s chord was not used for limb fitting since it had an 
inaccurate absolute time. 


be caused, for example, if the camcorder delayed to start 
the recording and, since this was an unattended pre-pointed 
station, this fact would not be noticed by the analysis of the 
video itself. 

The immersion recorded at Owings also seems to be 
shifted (delayed) with respect to the nearby chords (see Fig. 
[^. This chord, actually, has roughly the same length of Me- 
chanicsvllle’s, despite the fact that they are separated by 
about 100 km. Differently from Varina, though, this station 
had time stamps inserted in each frame of the video, which 
makes it more unlikely to justify an eventual timing issue. 
Another possible explanation for the apparent problem of 
the times of this chord is the determination of the ingress 
and egress instants in the light curve analysis, which could 
have been affected by noise. Finally, the delay during the 
Ingress could be caused by a relief feature in Ceres; we shall 
soon return to this hypothesis. 

In a second limb fitting, thus, we did not consider 
Brookline, Varina and Owings chords. The adjust of the 
hve parameters which define an ellipse to the 12 contacts 
then resulted in Xr,min ~ 1.27, indicating that the fitting 
is in good agreement with the observed data within the er- 
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Figure 5. The nine occultation light curves normalized and ver¬ 
tically shifted by a factor of 1.0 for better viewing, see Fig. 
for the explanation of the graph. The light curve of Brookline is 
shifted by -64 s as explained in the text. 


ror bars. This is the solution depicted in Fig. where we 
see that the chord length measured in Brookline is compat¬ 
ible to the model. The associated physical parameters are 
presented in Table as the nominal solution. 

For this event the polar aspect angle is = 90.7°, which 
makes the true oblateness equal to the apparent one, within 
the error bars. 

The position angle of this nominal solution is better 
constrained than the 2010 occultation. Actually, the uncer¬ 
tainty in the former is of 5° in contrast to 10° of the latter. 
This suggests that the pole-constrained solution obtained 
via the position of Ceres’ pole (Eq. would not be signifi¬ 
cantly different from the nominal one. 

This assumption was confirmed when we carried out the 
limb fitting with the constraint P = (25 ± 3)°, the position 
angle at the occultation which follows from Eq. The phys¬ 
ical parameters related to this pole-constrained solution are 
presented in the last column of Table and are essentially 
the same of the nominal solution. 

In the hypothesis that the delay observed in the immer¬ 
sion at Owings could be associated to a limb topography fea¬ 
ture, the recorded contact would then correspond to a neg¬ 
ative elevation of 31 ± 4 km with respect to the best-fitting 
ellipse. Theoretical models, however, predict that reliefs in 


Ceres should not be higher than about 10-20 km (John- 
[son fc McGetchin|1973| l, while published observational data 
sets the bound of 18 km (Carry et al.|200'8l. More recently. 


Ceres. 25 October 2013 



500 0 -500 

f,, (km) 


Figure 6. The best elliptical fit for the occultation chords for the 
event of 2013 using timing from Table|^and the pole-constrained 
solution. The arrow indicates the direction of motion, blue lines 
are the observed chords, the red and green segments are the 
ingress, egress and mid-occultation error bars at l-cr level. The 
chords with green error bars were not used during the limb fit 
process. The chord of Brookline if shifted by -64s as explained in 
the text. 


images by the probe Dawn also reveal an even smoother sur¬ 
face. Therefore, the association of Owings first contact to a 
relief is improbable. 


4 ASTROMETRY FROM OCCULTATIONS 


The objects’ geocentric positions derived from stellar occul- 


tatlons are most valuable for improving their orbits (Des- 


mars et al. 20151. Usually the precision of the astrometric 


positions are limited by the accuracy of the occulted star po¬ 
sition, not by the limb fit. Ceres’ geocentric J2000 positions 
at the time of each occultation are displayed in Eqs.j^and 
The errors of the positions come from the errors of the 
star positions, taken from the catalogues and from the errors 
of the relative apparent distances between star and Ceres, 
derived from the limb fit (which are displayed as Aq cos 5, 
AS). 


2010 Aug 17 < 


Time = 22 : 40 : 00 
a = 17‘'18“29’10122 ± 0'.'027 
S = -27°26'38'.'617 ± 0"028 
Aa cos S = 0'.'003; AS = 0'.'007 


(5) 


2013 Oct 25 


Time = 09 : 45 : 00 
a = 11*'57“52:9154 ± 0'.'019 
S = -b09°07'49'.'865 ± 0'.'021 
Aq cos S = 0"002; AS = 0'.'007 


( 6 ) 
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5 DISCUSSION 


A quick glance at Table shows an overall agreement be¬ 
tween the physical parameters derived from both occulta- 
tions, specially in the equatorial diameter. The differences 
of the solutions occur basically on the size of their error bars, 
and can be justified by the particularities of each set of data, 
as discussed below. 

The 2010 event, for example, had only seven contacts; 
none the less, they were well distributed over Ceres’ disc 
(see Fig.[^ acting as a constraint to its shape. On the other 
hand, the 2013 event had five more exploitable contacts, 
but they were concentrated In certain regions of the body. 
In particular, the absence of chords close to Ceres’ south 
pole made its oblateness less well determined here than in 
the 2010 event. 

However, even our best measurement for the oblateness, 
e = 0.08 ±0.03, has still a higher uncertainty with regard to 
other figures published in the literature, as Table shows. 
A larger number of uniformly spaced chords would be nec¬ 
essary to offer a best constraint to the oblateness. 

The few chords of the 2010 occultation could themselves 
only constrain the position angle of the object to a uncer¬ 
tainty of 10°. This uncertainty was reduced by a factor of 
two in the 2013 event, approaching -and verifying- the re¬ 


sult predicted by the work of Drummond et al. (20141. As 
was shown, using the coordinates of Ceres’ polar axis to 
limit the position angle was not an efficient procedure in 
the 2013 event, in the sense that it did not result in sig- 
nihcant changes in the parameters obtained in the nominal 
solution (see Table [^. 

On the other hand, constraining the position angle on 
the 2010 occultation was proved to reduce the error bars of 
the other parameters (disregarding oblateness). Moreover, 
this procedure resulted in an excellent agreement between 
the equatorial radius hgures of both events. 

The 2013 occultation, therefore, offers not only an in¬ 
dependent verification of the figures resulted from the 2010 
event, but also validates the procedure carried out there 
which led to the best-constrained parameters in this work. 

Comparison of Ceres’ equatorial diameter as measured 
by different techniques is carried out in Table We note 
an overall agreement between our result to those obtained 
via direct imaging by the Hubble Space Telescope (HST) 


(Thomas et al. 20051, the Keck Observatory and the ESO 
VLT ( Drummond et ^[2014 1. The smaller hgure reported 
by [Carry et al.| ( |2008| ) may be justified by the fact that in 
this study the effects of limb-darking were not taken into 


account, as pointed out by Drummond et al. (20141. 


As mentioned in the Introduction, the 1984 event (Mil 


lis et al. 19871 is the only occultation data to which we can 


compare our result. The measured diameters do not agree, 
ours being larger by 2-sigma. It is difficult to state for sure 
the reasons for this divergence. One way to clarify the issue 
is to redetermine the immersion and emersion instants from 
the original light curves using the same methodology pre¬ 
sented in our work. Note that in this occultation a variation 
of 0.5 s in the contact times correspond to almost 7 km on 
Ceres’ limb, which is on the same order of the residuals of 
their elliptical fit. However, we could not do that, as we had 
no access to the original 1984 light curve data. Moreover, 
we have no description of how the immersion and emersion 


Table 4. Ceres’ equatorial diameter and oblateness 


Eq. diameter (km) 

Oblateness 

Method 

Ref. 

972 ±6 

0.08 ±0.03 

Occultation 

1 

967 ± 10 

0.078 ±0.015 

Keck+VTL 

2 

959 ±5 

0.074 ± 0.007 

Keck 

3 

975 ±4 

0.067 ± 0.005 

HST 

4 

959 ±5 

0.05 ±0.01 

Occultation 

5 


References. 1: Pre sent work. 2: [Drummond et al.| (|2014[l. 3 : 
[Carry et al.| ( |200^ . 4: |Thomas et al.| ( |2005| (. 5: |Millis et al.| ( |1987| |. 


instants of the chords were actually determined from these 
light curves in Millis et al. (19871. 

NASA’s Dawn mission shall bring light to these ques¬ 
tions, which will be important not only for the knowledge 
of Ceres itself, but also for all the techniques used so far to 
study the physical properties of small Solar System objects, 
such as the stellar occultations. 
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